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^ '. Abstract 

^^ I The annihilation cross section of thermal relic dark matter determines both its relic density 

QQ ' and indirect detection signals. We determine how large indirect signals may be in scenarios with 

C^^ ■ Sommerfeld-enhanced annihilation, subject to the constraint that the dark matter has the correct 

.^ . relic density. This work refines our previous analysis through detailed treatments of resonant 

^y-T ! Sommerfeld enhancement and the effect of Sommerfeld enhancement on freeze out. Sommerfeld 

^D I enhancements raise many interesting issues in the freeze out calculation, and we find that the 

^"■^ ' cutoff of resonant enhancement, the equilibration of force carriers, the temperature of kinetic 

decoupling, and the efficiency of self-interactions for preserving thermal velocity distributions all 
play a role. These effects may have striking consequences; for example, for resonantly-enhanced 
Sommerfeld annihilation, dark matter freezes out but may then chemically recouple, implying 
C^ I highly suppressed indirect signals, in contrast to naive expectations. In the minimal scenario 

with standard astrophysical assumptions, and tuning all parameters to maximize the signal, we 
find that, for force-carrier mass m^ = 250 MeV and dark matter masses mx = 0.1, 0.3, and 1 
TeV, the maximal Sommerfeld enhancement factors are SgQ = 7, 30, and 90, respectively. Such 
boosts are too small to explain both the PAMELA and Fermi excesses. Non-minimal models may 
require smaller boosts, but the bounds on 5efr could also be more stringent, and dedicated freeze 
out analyses are required. For concreteness, we focus on 4/i final states, but we also discuss 4e 
and other modes, deviations from standard astrophysical assumptions and non-minimal particle 
physics models, and we outline the steps required to determine if such considerations may lead to 
a self-consistent explanation of the PAMELA or Fermi excesses. 

PACS numbers: 95.35.+d, 95.85.Ry 
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I. INTRODUCTION 

Dark matter may be composed of thermal relics with mass near the weak scale m^eak ~ 
100 GeV — 1 TeV. Such candidates are produced in the hot early Universe and then freeze 
out when the Universe cools and expands. Their annihilation cross section therefore deter- 
mines both the relic density and the rate of annihilation today. The requirement that the 
candidate be much or all of the dark matter therefore constrains its annihilation rate now, 
with important implications for indirect searches for dark matter. 

Although annihilation in the early Universe and now is determined by the same dynamics, 
the kinematics are vastly different: at freeze out, thermal relics have relative velocity frci ~ 
0.3, whereas today, v^ei ~ 10~^. As a result, the numerical values of the annihilation cross 
sections may differ significantly. For example, neutralino annihilation is in many cases 
dominated by P-wave processes [jj, and so is suppressed at low velocities. In this work, we 
consider Sommerfeld-enhanced cross sections, which have the opposite behavior: they are 
enhanced at low velocities and therefore boost present-day annihilation signals. We refine 
a previous study [2|] and examine how large these boosts may be, subject to the constraint 
that the thermal relic has the correct thermal relic density to be dark matter. 

Along with its general implications for future indirect searches, this study also has direct 
implications for the interpretation of current data. Following earlier excesses reported by the 
HEAT Collaboration [3l,3, the PAMELA Q, ATIC @, and Fermi Q Collaborations have 
reported excesses of cosmic positrons over an estimate of expected background M . Dark 
matter annihilation with Sommerfeld enhancement |9| of dark matter annihilatio n |lO| . gen- 



eralized to massive force carriers [11H15|, has been proposed as an explanation |16Hl9|. In 
its simplest form, this scenario assumes that dark matter is composed of a single particle 
species X, which interacts with light force carriers (p with fine structure constant ax and 
rrirj, <C mx ~ m^eak- This new interaction modifies dark matter annihilation and scatter- 
ing properties. The resulting annihilation cross section multiplied by relative velocity is 
icranVTe\)oS, where (cTan'^reOo IS its trec-lcvel value, and S is Sommerfeld's original enhance- 
ment factor [9|] 

1 _ Q-Tvax/v f ' 

generalized to massive 0, where v = frei/2 is the dark matter particle's velocity in the 
center-of-mass frame. In the proposed explanation, dark matter freezes out at early times 
when the Sommerfeld effect is negligible with (cTan'^reOo ~ 3 x 10~^^ cm^/s, leading to the 
correct relic density Q-ouh"^ — 0.114. At present, however, when v^ei is much smaller, the 
Sommerfeld effect becomes important, and, for example, for mx ~ 2 TeV and an assumed 
enhancement factor of 5 ~ 1000, such annihilations are sufficient to explain the positron 



excesses l2d, 21 



In Ref. [2(], we showed that in straightforward models, such large Sommerfeld enhance- 
ments cannot be self-consistently realized. The problem is simple to state: large Sommerfeld 
enhancement requires strong interactions, and strongly-interacting particles annihilate too 
efficiently in the early Universe to be all of the dark matter. More quantitatively, the re- 
quired new force carrier interaction necessarily induces an annihilation process XX — )■ 00, 

with cross section 

TTa\ , , 

(o-an^reljO — ■ (2) 

mx 
Conservatively neglecting other annihilation processes, for a typical weak-scale mass m,x, 



the thermal rehc density imphes a typical weak coupling ax ~ 0.05 \mxl'^ TeV], as dic- 
tated by the WIMP miracle. The resulting Sommerfeld enhancement now is ~ irax/v ~ 
100 [mx/2 TeV], an order of magnitude too small to explain the PAMELA and Fermi signals. 
Alternatively, to achieve S ~ 1000, one requires ax ~ 1, which implies Vtx ~ 0.001, two 
orders of magnitude smaller than the value assumed in deriving the requirement S ~ 1000. 
In this work, we refine our previous analysis in several ways. In Ref. |2| we approximated 
the Sommerfeld enhancement S by its value at m^ = 0, given by 5*° in Eq. ([T]). For 
massive (/), the Sommerfeld enhancement cuts off at a value proportional to mx/Tn^ and 
also exhibits resonant structure. Here we use numerical results (and a highly accurate 
analytic approximation to the numerical results 22h25|) for S. 



In addition, we refine our previous work to include the effect of Sommerfeld enhancement 
on freeze out. One might expect this effect to be negligible, because freeze out is typically 
thought to occur at f ~ 0.3, when the Sommerfeld effect is insignificant. However, anni- 
hilation continues to much later times, and Sommerfeld enhancement has an impact when 
the dark matter cools. Sommerfeld effects on freeze out have been considered previously in 



Ref. [26| for the case of Wino dark matter and in Refs. |27H30| for hidden sector dark mat- 
ter. In the analyses published after the PAMELA and Fermi excesses were reported |28l-l30|. 
((Tan'Wrci)o was taken as a free parameter, and the fact that it depends on the same parameters 
that determine 5" was not used. Here we make essential use of the observation that there is 
a irreducible contribution to (o"an'i^rei)o of the form of Eq. ([2]), and so constraints on the relic 
density bound the fundamental parameters that also determine S. 

As we will see, the combination of these two refinements highlights a number of effects 
that may typically be ignored, but now require exploration. For example, the relic density is 
in principle affected by the cutoff of resonant Sommerfeld enhancement, the production and 
decay rates of force carrier particles, the temperature of kinetic decoupling, and the efficiency 
of self-scattering for preserving the dark matter's thermal velocity distribution. Very close 
to resonances, we also find the intriguing possibility of chemical recoupling, in which dark 
matter freezes out and then melts back in, with annihilations becoming important again 
at late times. This implies that exact resonances suppress, rather than enhance, indirect 
signals. 

We examine the quantitative impact of all of these effects, and present results for the 
maximal Sommerfeld enhancement of indirect signals. In the most optimistic case, that 
is, tuning all parameters to maximize the signal, we find that, for m^ = 250 MeV and 
nix = 0.1, 0.3, and 1 TeV, the largest possible effective Sommerfeld enhancements are 
Ses = 7, 30, and 90, respectively. This refined analysis therefore strengthens our previous 
results: Sommerfeld enhancements may be significant, but they fall short of explaining the 
current PAMELA and Fermi cosmic ray excesses. We then discuss various astrophysical 
effects that may reduce the discrepancy between the maximal enhancements derived here 
and those required to explain the data, and critically examine several non-minimal models 
and the issues that must be addressed to determine if more complicated particle physics 
models may provide viable explanations. 

This paper is organized as follows. In Sec. [ITl we describe our underlying model as- 
sumptions and our treatment of Sommerfeld enhancements with resonances. In Sec. IIIII we 
analyze the effect of Sommerfeld enhancement on freeze out. In Sec. [IV]we assemble these 
pieces and present the results for the maximal Sommerfeld enhancement achievable now. In 
Sec. lYl we compare these to the current data, and discuss non-standard astrophysical effects 
and non-minimal particle physics models. We present our conclusions in Sec. IVI[ 



II. SOMMERFELD ENHANCEMENT WITH RESONANCES 

We consider a simple model with a hidden sector dark matter particle X, which cou- 
ples to a light force carrier cj) with coupling y/Anax- The annihilation cross section is then 
(cTan'^rcOo'S', where (cranfrei)o IS the tree-level cross section and S is the Sommerfeld enhance- 
ment. 

To maximize the Sommerfeld enhancement, we assume that (cranfrei)o is dominated by 
S'-wave processes, and so is unsuppressed at low velocities. We also consider only the 
"irreducible" annihilation channel XX — )■ dch, and take the tree-level cross section 



') 



Tiax 



O"anfrcl)o = o" " (^) 



ni 



x 



This may be modified by 0{1) pre- factors, depending, for example, on whether X is a 
Majorana or Dirac fermion, and whether is a scalar or a gauge boson. In addition, 
even in simple models, (cranfrei)o will typically receive additional contributions from other 
annihilation channels; to maximize the Sommerfeld effect on indirect search signals, we 
neglect these other channels here, but discuss their impact in Sec. |Vl 

To determine the enhancement factor S, we numerically solve the differential equation 

1 d\ ^ax^_^^r^^_^^^2^^ (4) 



nix dr"^ ^ 

with the boundary conditions x'{^) = irrixvxir) and x(r) = e*"^-^'"' when r — t- oo. The 
Sommerfeld enhancement factor is given by 

5 = M^. (5) 

The Sommerfeld enhancement may also be obtained by approximating the Yukawa po- 
tential by the Hulthen potential, for which an analytic solution is possible 1^24,]. The resulting 



analytic approximation to the Sommerfeld enhancement is [24, [25 
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S = / ^^ "^ ^=^ , (6) 
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where e„ = v/ax and e^ = m^/{axmx)- The analytic expression of Eq. ([6]) is compared 
to the numerical solution in Fig. [H We see that the analytic result is an excellent approxi- 
mation, typically reproducing the numerical results to within fractional differences of 10%, 
and accurately reproducing the resonant behavior. Given Eq. ([6]), we see that for e^ ^ e,/,, 
these resonances are at 

^<^- — S"^~ ' ^ = 1,2,3... (7) 

At these resonances, with m^ determined by Eq. (I7l), the Sommerfeld enhancement for low 
V is 

S:.^-^. (8) 

S is therefore enhanced by w~^ at resonances, as opposed to v~^ away from resonances. 
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FIG. 1: The Sommerfeld enhancement factor 5 as a function of e^ = m(j)/{axfnx) for the constant 
values of e^ = v/ax indicated. The sohd red curves are the analytic approximation of Eq. ([6]), and 
the dashed blue curves are numerical results. 



Note, however, that both the numerical and analytic results become infinite on resonance. 
This is unphysical, a result of the fact that the quantum mechanical treatments do not 
include the effect of bound state deca^ In fact, the finite hfetime of bound states implies 
that 5" saturates at f ~ a|^(m<^/mx) [l3|, which, given the condition for a strong resonance. 



Eq. ([7]) with small n, is v 
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In this study, we use the analytic form for S given in 



Eq. ([6]) and include the saturation by making the substitution f — )■ f + a 
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III. THERMAL FREEZE OUT WITH SOMMERFELD ENHANCEMENT 



A. General Formalism 



Given values of mx-, Tri<f>i and ax-, we must determine the relic density Vtx- The formalism 



of thermal freeze out is well developed |31|, |32| . Here we review this formalism in sufficient 



generality to accommodate novel effects resulting from Sommerfeld enhancements. 

The evolution of the abundance of a thermal relic X is governed by the Boltzmann 
equation 

^ + ?,Hnx = -(cTant^rcl) (^X " ^X ^) ) (9) 



dt 



where nx is the number density of the dark matter particles, rfx is its value in equilibrium, 
H is the Hubble constant, and (cTan'^rei) is the annihilation cross section multiplied by the 
relative velocity, averaged over the dark matter velocity distribution. Changing variables 
from t — >■ a; = rrix/T and nx -^ Y = nx/s, where s is the entropy density, the Boltzmann 
equation becomes 



dY 
dx 



vr {.9*slJ 

45 x^ 



■{cr^nVrel)(Y^-Y'^^), 



(10) 



where mpi = Gj^ ^ 1.2 x 10 GeV is the Planck mass, and (7*3 and g^ are the effective 
relativistic degrees of freedom for entropy and energy density, respectively. In this section 



we will assume that the particles X annihilates to are in thermal equilibrium at the time 
of freeze out. We will revisit this issue in Sec. IIIIBt where we show that this requirement 
leads to significant constraints if the dominant annihilation is to the dark force carriers. 

To evaluate the relic density, we first evolve it to the time of chemical decoupling, when 
the annihilation rate Fan = nx (can^^rei) is approximately the expansion rate H. As discussed 
in Sec. HIl we assume the annihilation cross section has the form (Tani'rei = (can'^reOo'S', where 
((Tan'^reOo IS the S'-wavc, tree-lcvcl cross section, and S is the Sommerfeld enhancement. Up 
to the time of freeze out, the Sommerfeld effect is insignificant, and so the standard results 
for S-wave annihilators apply: defining freeze out by 

Y{xj) = (1 + c) Y^^ixf) , (11) 

where c is a constant, freeze out occurs at 

Xf ^ In^- -ln(lnO 
^ = 0.038c(2 + c)mpimx (fi'/y^) (o-an^^roi)o , (12) 

where g is the number of degrees of freedom of the dark matter particle; we take g = 2. To 
match numerical results, c ~ 1; we choose c = 1 and have checked that varying c from 1/2 
to 2 has no appreciable effect. To further test this approximation we have also numerically 
evolved Eq. (ITOl) from x = 20 to a; = 100 (when the equilibrium abundance is effectively 
zero) and then used Eq. f fT2|) to evolve into the late Universe. This numerical calculation 
typically results in 10% smaller relic abundance. Given the excellent agreement, we set c = 1 
and with this choice, Eqs. fITT]) and flT2|) determine Y{xf), the abundance at freeze out. 

After freeze out, Y^^ quickly becomes insignificant. Neglecting it, we may solve Eq. f lTO]) 
to find 

1 1 nr Pl^d (5(^^/y^)(cranWrel) , 

+ \ h^iTT'Pi'mx — o dx 
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VT T' i9*s/^/g^){(^llnVTel) 



+ ^^mp, mx l^^ -, dx , (13) 

where Xkd is the value of x at kinetic decoupling and Xg is its value when annihilations 
become insignificant and we may stop the integration. We have broken the integral into two 
to emphasize that there are two eras: before and after the temperature of kinetic decoupling 
Tkd- Before T^d? the dark matter's velocity distribution is thermal with temperature Tx = T. 
After Tkd, the dark matter's velocity distribution initially remains thermal, but Tx drops as 
a~^, while T drops as a~^, where a is the scale factor, and so Tx = T^/T\^d- Eventually, the 
dark matter's velocity distribution need not even be thermal. We discuss the value of Tkd 
and the issue of non-thermal velocity distributions in Sees. IIIICI and IIIIDt respectively. 

When the dark matter distribution is thermal with temperature Tx, the thermally- 
averaged cross section in the non-relativistic limit is 

(o-an^rcl) = / f (Vl) f{v2)a^nVre\d^Vi(fv2 

2 1 — "'''^' 

:TV^^ye '^ as^nVreldVrel , (14) 
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where f{v) is the dark matter's velocity distribution, v^ei = \vi — V2\, and 



v„ = J?^^ J^ , (15) 

V mx ^ xx 

the most probable velocity. For the Sommerfeld-enhanced annihilation cross section, this 
becomes 



.3/2 ,^ 



where 



(cTan^^rel) = TT^ / (o-an^^reO^reie ''"'"'""'^^^dVrel = (<7an^^rcl)o5'(xx) , (16) 

2^yr^ Jo 

S{xx) = TTT / Sv^.e-^-'^'^^/'dv.ei (17) 

2y/TT Jo 

is the Sommerfeld enhancement averaged over a thermal distribution with temperature Tx = 

mx/xx- 

B. Equilibration of Force Carriers 

In the relic density calculation, we have implicitly assumed that, when the X particles 
chemically decouple, the force carriers (p are in thermal equilibrium with the massless par- 
ticles of the standard model. The force carriers are expected to interact with the standard 
model, as their decays to standard model particles typically provide the indirect signals. In 
a simple example, if is a U(l) gauge boson, it may mix with the standard model photon 
through kinetic mixing terms ~ eF^^F"^^^ . After diagonalizing the 0-photon system, stan- 
dard model particles with charge Q have hidden charge eQ J33|, and so the particles decay 
through — )• //, where / = e, /i, m, d,s, . . .. The largest kinetic mixing parameter allowed 
by current particle physics constraints is e ~ 10~^ 3^, |35| . 



The existence of interactions with the visible sector does not, however, guarantee that 
they are efficient enough to bring the particles in thermal equilibrium. Here we determine 
sufficient conditions for the kinetic mixing example to guarantee the equilibration of 
particles. The leading number-changing interactions between the particles and the visible 
sector are decays — )■ // and inverse decays // -^ 0. Other number-changing processes, 
such as // — )■ 07, qq — )■ (j)g, and f'j — )■ /0, are parametrically suppressed compared to 
these and subdominant. For temperatures T ^ m^, the decay and inverse decay processes 
balance (since /, / are in chemical equilibrium), and we may simply compare the decay rate 
to the expansion rate |36| . 

For T ^ m^, the thermally-averaged decay rate for — )■ // in the lab frame is (F,^) ~ 
(e^/3) J2f{Q}Nl)aEM^TT'<t>{^TT'(t>/T), where we have averaged over three polarizations, Qf and 
N^ are the standard model charge and number of colors for fermion /, the sum is over all 
fermions with mass rrif < m^/2, and the last factor oim^/T is from time dilation. Note that 
the decay process is inefficient at early times and high temperatures, but as T drops, (F^) 
increases and H decreases, and so the particles may come into thermal equilibrium with 
the the standard model. The expansion rate is H c::^ 1.66 gj T^/rripi; at freezeout, g^ ^ 100. 
Requiring (F,^) > if at Ty ~ mx/25 yields e ^ 3 x 10~^ for m,^ = GeV and mx = TeV. 

If the above constraint is not satisfied, then the hidden sector may be at a different 
temperature from the standard model, modifying the dark matter freezeout. Such scenarios 



have been considered in Ref. J37|], and the difference in temperatures has a significant impact 



on the dark matter relic density. These considerations may be important for some models, 
but here we assume that Eq. (fTSi) is satisfied and m^ <^Tf '2:i mxl'2.^, so that the particles 
are in thermal contact with the standard model at freeze out and their number density is 
close to the thermal prediction. 

We note here that the calculation above does not guarantee chemical equilibrium. For 
the particles to produce X through cfxp — )■ XX, we need the tail of the distribution to 
be populated thermally. To check this, we approximate the decay rate for with energies 
of order mx as (e^/3) X]/(<5/^/)ttEM^0('TT'0/^x)- If this decay rate is larger than the 
expansion rate, then the inverse process will also be in equilibrium and hence thermally 
populate the tail of the distribution. Comparing this to the expansion rate at Tf ~ mxj'^^ 
we obtain that 



e > 1.5x 10" 
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(18) 



In setting the bound above, we have not considered scattering interactions with standard 
model fermions (and the resulting changes in kinetic energy for 0) because these processes 
are proportional to e" and hence considerably slower that the inverse decay process. 

To summarize, for the particles to be in thermal equilibrium with the standard model 
at freezeout requires that Eq. (ITS]) is satisfied and m^ <t^Tf c^ mx/25. We will assume these 
conditions hold in our analysis, but we note that they do not necessarily hold. In particular, 
if m^i, ~ Tf, the number of particles at freezeout is reduced, which reduces the X relic 
density, strengthening the bounds on Ses we determine below. 



C. Kinetic Decoupling of Dark Matter 



After freeze out, dark matter particles remain kinematically coupled to the thermal bath 

through elastic scattering. Kinetic decoupling occurs later, when the momentum transfer 

rate drops below the Hubble expansion rate. We define the kinetic decoupling temperature 

Tkd by rfc(Tkd) = ii/^(Tkd), where r^ is the momentum transfer rate. It may be approximated 

as 

T 

Tk ~ UriaelVrel) , (19) 

mx 

where rir is the number density of the relativistic species in the thermal bath, and (cieifrei) 
is the thermally-averaged cross section for elastic scattering between dark matter particles 
and the relativistic particles in the thermal bath. 

We first consider the model-independent elastic scattering off the hidden sector thermal 
bath of particles (38|. For T < m^, n^ is Boltzmann suppressed, and so this process will 
not be able to maintain kinetic equilibrium. For T ^ m,^, the thermally- averaged X0 — )■ X0 
cross section is 



(o-cl^^rcl) 



a 
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the Thomson scattering cross section. Given the number density n^ 
transfer rate is F^ ~ a'j^T'^/m,^, and this is equal to iJ ~ T^/mpi at 



(20) 
T^, the momentum 
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(21) 



where we have normahzed the expression to typical parameters that give the correct rehc 
density in the presence of Sommerfeld enhancement. Halo shape constraints require m^ ^ 
30 MeV [2:], and so for all physically viable and relevant parameters, 



T, 



kd 



m^ . 



(22) 



This is the temperature of kinetic decoupling from the hidden thermal bath. It is quite 
model-independent, as it assumes only the X-cp interactions required in all Sommerfeld 
enhancement scenarios, and so it provides a maximal value of T^d in Sommerfeld scenarios. 
Dark matter particles may also be kept in kinetic equilibrium through scattering off the 
visible sector's thermal bath. This is more model dependent, but as discussed in Sec. lIIIBl if 
is a U(l) gauge boson, it may mix with the standard model photon. After diagonalizing the 
0-photon system, standard model particles with charge Q have hidden charge eQ, inducing 
new energy transfer processes. The most efficient process is Xe — )■ Xe scattering through 
t-channel exchange. For T < rUe, Boltzmann suppression of the electron number density 
makes the interaction inefficient. At temperatures m^ ^ T ^ rUe, the corresponding cross 



section is 39 
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(23) 



The momentum transfer rate F^ ~ e^Q;EM«x^^/("^l^™x) is equal to the expansion rate at 
temperature 
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nl/4 



(24) 
e'aEMax'o^pi. 

and so the resulting temperature of kinetic decoupling from the visible sector's thermal bath 
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(25) 



where we have again normalized ax and mx to typical freeze out parameters, and addition- 
ally normalized m<^ to its smallest possible value and e to its largest allowed value. When 
e is near its maximal value, T^^{t) < Tj^, and so interactions with the visible thermal bath 
delay kinetic decoupling. We denote the lowest possible kinetic decoupling temperature 



T, 



kd 



^kdle 



(e = 10-3) 



(26) 



and will explore the dependence of our results on the temperature of kinetic decoupling by 
varying it between its maximal value T^^ and its minimal value T^^. 

By crossing symmetry, the visible sector scattering interaction also implies an annihilation 
process XX — )■ e^e~ . In principle, this also enters the thermal relic density calculation. 
The cross section is (cTan'^rei) ~ e^ctEMCtx/""^!:, however, and the e^ suppression makes this 
subdominant for all relevant cases. 



D. Velocity Distribution after Kinetic Decoupling 

Before kinetic decoupling, dark matter particles have the same temperature as the thermal 
bath and a Maxwell-Boltzmann phase space distribution. Usually one assumes that the 



9 



phase space distribution remains Maxwell-Boltzmann after kinetic decoupling. However, 
this is not necessarily true in scenarios with Somnierfeld-enhanced annihilation, because 
slow particles annihilate with larger cross sections. This preferentially depletes the low 
velocity population and may distort the phase space distribution. In this case, the standard 
formulae used to compute the relic density of dark matter are not valid and one needs to 
explicitly consider the effect of annihilations using the full phase space distribution. 

As discussed in Refs. [2|, |39|-l42|, however, the field mediates self-scattering between 
dark matter particles. If the self-scattering rate is higher than the Hubble expansion rate, 
the momentum exchanged in these self-interactions with be sufficient to maintain thermal 
equilibrium. The momentum transfer of particles interacting through Yukawa potentials 
has been studied in Refs. |43l . |44| . Although the authors of these studies were interested 
in slow and highly charged particles moving in plasmas with screened Coulomb potentials, 
they approximated these potentials by Yukawa potentials, and so their results are exactly 
applicable in the current particle physics context. The numerical results of these studies for 
the momentum transfer cross section are accurately reproduced by |43l . |44 
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where /3 = 2ax^(/>/(^x'yrei)- 

Given these results, we may calculate the dark matter self-scattering rate and compare 
it to the Hubble expansion rate. The self-scattering rate is given by [2| 



f{Vl)f{v2) rix (fXTWrel) ^d^Vid^ 



(28) 



where nx is the dark matter density. This sets the time scale to change velocities by 0{1). 
To find the rate to change the kinetic energy by C(l), one would divide the above rate by a 
factor of 3. However, such details are not important for the present calculation. H and Tg 
are presented in Fig. [2] for two representative cases. 

We may define the self-scattering decoupling temperature T^t by rs(Tnt) = H{T^t)] after 
Tnt, the dark matter velocity distribution may become non-thermal. To derive an approxi- 
mate expression for Tnt, note that to have the right relic abundance, the co-moving number 
density of dark matter during the self-scattering epoch cannot be smaller than the co-moving 
number density at present. The dark matter number density nx may therefore be taken to 
be 



nx ^ s 



(nx\ 



^DMPc 



s X 3.9 X 10" 



-13 



TeV 



(29) 



V s /o rnxSQ I nix 

During the self-scattering decoupling epoch, typically f3 > 1000, and so ax ~ i^/fn% where 
K ~ 7rln^/3 ~ 600 — 1000 for the parameters of interest. With this approximate expression 
for ax, we find 
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FIG. 2: The self-scattering rate Tg (solid red) and the Hubble rate H (dashed blue) as functions 
of j; = mx/T for the values of mx, m^, ax, and T^d indicated. 

This agrees well with our numerical results, and we use this approximate form with k = 800 
in deriving our results below. 

From both the numerical and analytical results, we see that the self-scattering rate effi- 
ciently preserves the thermal distribution to temperatures T^t ~ 10 keV, when dark matter 
velocities are v^ei ~ T/(Tkd?7ij5f)^^^ ~ 10~^. This is a low velocity, and in most cases, we will 
find that it is sufficiently low that the impact of non-thermality after T^t has a negligible 
impact on the thermal relic density. Note, however, that Tnt is not low enough to ensure a 
thermal distribution down to Vj-ei ~ a^, where the Sommerfeld resonances cut off, and so 
very close to resonances, the non-thermality will have an effect. 



IV. RESULTS 



A. Maximal Sommerfeld Enhancements 



In this section, we present results for the maximal Sommerfeld enhancement. We have 
defined several Sommerfeld factors, including S^, the original Sommerfeld factor for massless 
(j) given in Eq. ([T]), and S, the Sommerfeld factor generalized to massive 0, which includes 



resonances and is given in Eqs. (jS]) and 
enhancement factor 



Here we also define the effective Sommerfeld 
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Xnow = 2/^0, and A^ = eTi{z/y/2) — (2/7r)^/2^g~2 /2 -^jth z = fmax/'^o- •S'eff is the experimen- 
tally relevant parameter, as it is the factor by which indirect fluxes are enhanced relative to 
the case without Sommerfeld enhancement. Unlike the early Universe case, in the halo we 
have to cut off the velocity integral at some maximum relative speed fmax? which is related 
to the escape speed from the local neighborhood. This maximum speed is a function of the 
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angle between the dark matter particles in the lab frame. However, we have checked that 
setting fmax to be equal to the escape speed fesc and integrating as in Eq. (132|) makes a 
difference of less than 10%. The escape speed at about 10 kpc from the center of the halo 
is estimated to be about 500 km/s to 550 km/s |45[. We use fesc = 525 km/s to derive our 
results, noting that 10% variations in fgsc have a much smaller effect on See- We have as- 
sumed that the velocity dispersion tensor is isotropic with the 1-D velocity dispersion given 
by Wo/\/2. We use vq = 210 km/s consistent with determinations of the circular velocity 
(e.g., |45|). Inferring the dispersion from the circular velocity curve requires knowledge of 
the dark matter density profile. The value we use for vq is consistent with a Navarro- Frenk- 
White density profile with a scale radius of about 20 kpc. Here again, we note that there are 
uncertainties in the value of fo at the 10% level. This translates directely into an uncertainty 
in Ses of about 10% if S' oc 1/v and 20% around resonance where S oc 1/f^. 
For comparison purposes, we also define 

^^^^ " 3 X 10-26 cm3/s ' ^'^'^' 

where 

^°ow = ^^ I S%l,e-^-^^^^"dv,e, (34) 

is the effective Sommerfeld enhancement without resonances, but with the Sommerfeld effect 
on freeze out included. 

To calculate the largest possible S'eg for a given irix and m^, we make a number of 
optimistic (S'eg-maximizing) assumptions: 

• We fix Tkd = Tj^j. This delays kinetic decoupling as much as possible, keeps the 
dark matter as hot as possible, and so reduces the Sommerfeld effect on freeze out, 
maximizing S^q. 

• We fix Xs in Eq. (IT^ by stopping the dark matter evolution at T^t] that is, we neglect 
all annihilations that occur after the dark matter distribution becomes non-thermal. 
This is certainly optimistic, as the distribution will remain thermal for some time 
and dark matter annihilations will continue, but this again minimizes the Sommerfeld 
effect on freeze out, maximizing Ses- 

• We require Qxh"^ = 0.114. This might appear to be too restrictive; after all, there is 
no requirement that the observed signals arise from a particle that makes up all of 
the dark matter. However, if a flux arises from Sommerfeld-enhanced annihilation, 
it scales as n^(aanWrei)'S' ~ ax^, because its number density scales as nx ~ ^x ^ 
(o'anVrei)'^ ~ «x^ and S ~ ax- The flux can therefore always be increased by lowering 
ax until Vtxh"^ is all of the dark matter, and so choosing Vlxh"^ in fact maximizes 
indirect signals. 

• We choose the maximal ax that yields Vtxh'^ = 0.114. Roughly speaking, Qxh"^ 
decreases as ax increases, and so typically, there is a unique choice of ax that yields 
the correct Qxh^- In some cases with strong resonances, however, there are three 
choices of ax that give the correct Qx , as shown in Fig. [3J When we do a coarse scan 
there is no guarantee that we always pick the solution with the largest value of Ses- 
Note, however, that the values of Scs only change by about 20-30% between the three 
allowed solutions, and such variations do not modify our conclusions. The values of 
ax that yield the correct relic density for given values of mx and m^ are given in 
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FIG. 3: The relic density Q.xh? (solid red) and S^s (dotted blue) as a function of the fine-structure 
constant ax for the fixed values of mx-, m^, and T^d indicated and the observed value of Q.xh? 
(dashed black). In most cases, there is a unique choice of ax that yields the correct Qx^"^ = 0.114 
(left). In the presence of strong resonances, however, there are cases where three different choices 
of ax all yield the correct Qx (right). In these cases, Ses varies by about 20-30% between the 
different solutions. 
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FIG. 4: The value of ax required to achieve a relic density of Qxh"^ = 0.114 as a function of the dark 
matter particle mass nix (solid red) for ?7i^ = 250 MeV (left) and m^j, = 1 GeV (right). We also 
plot the required ax (dotted blue) if Sommerfeld effects are neglected in the early Universe. The 
tree level cross section (without Sommerfeld enhancement) is ((Tan't'rei)o = ira'j^/mj^. Because ax 
varies over almost two orders of magnitude in this plot, the dips near resonance are not immediately 
apparent. 



Fig. m Note that close to a resonance, one needs extreme fine-tuning to avoid efficient 
annihilation in the early Universe and obtain the correct relic density. 

In Fig. we show the maximal values of S^s as a function of mx for m^ = 250 MeV and 
1 GeV. To understand the impact of resonances and the Sommerfeld effect on freeze out 
on these results, we also plot other Sommerfeld enhancements. 5", which includes neither 
resonances nor freeze out effects, was used in Ref. (2(]. We see that it is almost always 
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FIG. 5: The effective Sommerfeld enhancement factor S^s (sohd red) as a function of nix for 
m^ = 250 MeV (left) and 1 GeV (right) and the set of Sefj-maximizing assumptions hsted in 
Sec. IIV A[ Also shown for comparison are S^ (dotted blue), the Sommerfeld factor without reso- 
nances and neglecting the Sommerfeld effect on freeze out; S (dot-dashed blue), the Sommerfeld 
factor with resonances but neglecting the Sommerfeld effect on freeze out; and 5gg (dashed red), 
the Sommerfeld factor without resonances, but including the Sommerfeld effect on freeze out. 



overestimates the maximal Sommerfeld enhancement. This is because, as evident in Fig. |5l 
the effect of Sommerfeld enhanced annihilation on freeze out is a significant suppression. 
This may be understood as follows. If Sommerfeld enhancement reduces the thermal relic 
density by a factor (, the tree- level cross section ((Jan'^reOo must be reduced by ( to keep fix 
fixed. However, Scs oc a^ oc (o"anfrei)o ; cind so reducing fix by ( implies a reduction in S^s 
by a factor C^^'^. For example, for mx = 1 TeV, the maximal coupling is ax — 0.021, and 

Including the effect of freeze out here 
2, a significant reduction. 
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C ~ 1.5, consistent with the results of Refs. |29l . 
therefore reduces 5° to S^^ by a factor of (1.5)^/^ 

Adding resonances then produces oscillations about the S^ and 5'°^ contours, with peak 
positions given by Eq. ([7]). We see that, although resonances can significantly enhance the 
effective Sommerfeld factor, these enhancements are significant only for low values of n; for 
m^ < GeV, these he at low mx, and the effect becomes negligible for mx ^ 1 TeV. 

In Fig. m we have fixed m^ to typical values considered in the literature. In Fig. |6l we plot 
upper bounds on Ses as a function of m,^ for fixed values of m^x- We see that the resonances 
have little impact for m,^ ^ 1 GeV, but can produce enhancements by factors of 2 to 3 for 
larger m^. Large values of m,^ > 1 GeV have been considered disfavored, however, as they 
eliminate the kinematic suppression of anti-proton fluxes, which are typically considered to 
be consistent with astrophysical backgrounds. 

What is the effect on Ses of deviating from the S'efr- maximizing assumptions? The effect 
of the choice of Tkd is highly sensitive to how close one is to a resonance. Far from a 
resonance, there is little sensitivity to Tkd values between T^ and T^^. As one approaches a 
resonance, however, Ses may vary by ~ 10% or more. In particular, we found regions close 
to resonances for tti^ ~ GeV where Ses changed by about 30%. This behavior is shown 
in Fig. [71 If one is not very close to a resonance, the choice of Xg has a small impact. In 
contrast, as discussed above, the assumption of flxh"^ = 0.114 has a large effect; for other 
choices, the maximal value of S^g scales as {flxh'^Y^'^. Finally, as can be seen in Fig. |31 
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FIG. 6: The effective Sommerfeld enhancement factor S'cfj as a function oim^ for the mx indicated 
and the set of 5efj-niaximizing assumptions hsted in Sec. IIV A[ 
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FIG. 7: The effective Sommerfeld enhancement factor 5ofr for T^j = Ty^^ (sohd red) and Tkd = ?kd 
(dashed blue), for ra^ = 250 MeV (left) and 1 GeV (right). All other 5efr-maximizing assumptions 
have been made. Near resonances, 5eff is sensitive to the temperature of kinetic decoupling. 

different choices of ax also decrease 5'efr oc oCx by ~ 20 — 30%. 

For typical parameters that are not very close to a resonance, then, we expect that 
assumptions different from those listed above would produce a small decrease in S'eg when 
away from resonances. Note, however, that very near a resonance, all of these choices become 
very important, as we now discuss. 
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FIG. 8: The annihilation rate Fan = rix {cr anVrei) (sohd red) and the Hubble rate H (dashed 
blue) as functions of x = rrix/T for the values of mx, rn^, ax and T^d indicated. Without 
resonant Sommerfeld enhancement (left), Fan < H after freeze out; but -with resonant Soninierfeld 
enhancement (right), Fan may become comparable to H at late times, leading to chemical recoup ling 
and an ne-w era of annihilation. Also indicated on the plots are the resulting values of Qxh'^ that 
result from assuming annihilations up to a temperature of Tnt -when the dark matter distribution 
function is no longer able to maintain kinetic equilibrium, as -well as a lower temperature of 1 eV. 

B. Resonances and Chemical Recoupling 



Very near a resonance, the annihilation cross section becomes sensitive to details of the 
dark matter's evolution to lo-w velocities at late times. The effective Sommerfeld enhance- 
ment then becomes highly sensitive to details of T^d, the velocity distribution of dark matter 
after Tm, and the cutoff of resonant enhancements. As this is far from generic, we do not 
present details of these dependences, but we note that, contrary to naive expectations, Ses 
is not maximized by sitting exactly at resonance. In fact, exact resonances lead to extremely 
efficient annihilation in the early Universe and minimize current indirect signals. 

Resonant Sommerfeld enhancement's effect on freeze out can also be so large that it leads 
to the intriguing phenomenon of chemical recoupling. This is illustrated for two represen- 
tative cases in Fig. |8l Without resonant Sommerfeld enhancement, dark matter freezes out 
and remains frozen out, but with resonant enhancement, dark matter may melt back in at 
late times, or chemically recouple, leading to a second era of efficient annihilation. 

The phenomenon of chemical recoupling in the context of Sommerfeld enhancements may 
be understood as follows. The Hubble parameter scales as if ~ T^. The annihilation rate 
is Fan = nxicTanVrei)- The dark matter density scales as nx ~ T^, and (cTan'^rei) scales as 
vZl ^ Tv off resonance and v~J ~ Tx^ on resonance. Before Tkd, Tx = T, and so the 
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annihilation rate cannot grow relative to the expansion rate. However, after Tkd, Tx oc T^, 
and so Fan ~ T. At late times, then, the annihilation rate decrease slowly enough to become 
comparable to or greater than the expansion rate. The dark matter then melts back in, and 
there is a new era of annihilation. In these cases, the relic density is very sensitive to the 
temperature at which we stop including the annihilation process. As our default case, we 
have been conservative and stopped the annihilation process when the self-interactions of 
dark matter are no longer able to maintain kinetic equilibrium. Annihilations will proceed 
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FIG. 9: The maximal effective Sommerfeld enhancement factor Sqs compared to data for force 
carrier mass m,^ = 250 MeV. The cross-hatched region is excluded by requiring consistency with 
the thermal relic density and adopting all of the S'efj-maximizing assumptions listed in Sec. IIVAI 
The red and green shaded regions are 2a PAMELA- and Fermi-favored regions for the 4/i channel, 
and the best fit point is {mx,Ses) = (2.35 TeV, 1500) fiF 



beyond this temperature; however, we cannot assume a Maxwelhan distribution for the dark 
matter particle, because the annihilations preferentially deplete the low momentum tail. As 
an example, we show in Fig. |8] what happens when we allow the annihilation process to 
proceed down to 1 eV. The relic density is essentially negligible for the case where we have 
chemical recoup ling, whereas it is unchanged for the case away from resonance. 



V. COMPARISON TO PAMELA AND FERMI 



A. Maximal Enhancements and Best Fit Parameters 



In Figs. M and [TOl we compare the maximal Ses presented in Fig. |5] to the boost factors 
required to explain PAMELA and Fermi data. The PAMELA and Fermi regions, derived 
in Ref. J20|, assume an isothermal halo and a 250 MeV force carrier that decays with 100% 
branching ratio to /i"*"/^", leading to annihilations XX — t- fi^fi~^^fi~. The fits are insensitive 
to m0, provided m^ <^ mx- The best fit is for {nix, Ses) = (2.35 TeV, 1500), and the regions 
are 2a contours relative to the best fit parameters. Results for other halo profiles, other 
particle physics models, and other final states have also been considered [20|, |2l|, |46| . These 
studies find that the 4/i final state provides a better fit than all other considered final states, 
but see Sec. IV Dl for a discussion of this issue. 

Despite choosing T^d and all other parameters to maximize Ses, we find that the possible 
values of Ses fall short of explaining the PAMELA and Fermi excesses. For example, at 
TUx = 2.35 TeV the maximal Sommerfeld enhancement is ^ofr ~ 100, a factor of 15 below 
the best fit value. For lower mx, the data may be fit with smaller Scs, but the relic density 
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FIG. 10: The maximal effective Sommerfeld enhancement factor Scs compared to data for m^ = 
1 GeV. The cross-hatched and shaded regions are as in Fig. [9l The PAMELA- and Fermi- favored 
regions are for m,^ = 250 MeV, but do not change appreciably for m,^ = 1 GeV [211] in the case 
that only the 4/x channel is included. 



bounds are also stronger; for example, for m^ = 250 MeV and rrix 
the maximal values of Ses are 7, 30, and 90, respectively. 



0.1, 0.3, and 1 TeV, 



B. Astrophysical Uncertainties 

Our analysis so far has assumed standard astrophysics and astroparticle physics. To 
enable Sommerfeld enhancements to explain PAMELA and Fermi, one may consider uncer- 
tainties in these assumptions. 

The annihilation signal is sensitive to the square of the dark matter density px in the 
local neighborhood. The best fit regions given in Ref. 20l| assume the conventional value 
Px = 0.3 GeV/cm^, which has traditionally been considered to be uncertain up to a factor of 
2 47| . More recently, studies have tried to refine the calculation of the uncertainty in the local 
density of dark matter. One study found px = 0.389±0.025 GeV/cm^ J48|, using an ensemble 
of kinematic tracers to constrain a standard galactic model with a spherical halo, stellar and 
gas disk and stellar bulge. Another recent study |49| finds px = 0.2 — 0.4 GeV/cm^. A 
third study Taylor-expanded the local rotation velocity curve to first order and used local 
kinematic tracers and an estimate of the baryonic contribution to the rotation curve to 
find Px = 0.430 ± 0.113 ± 0.096 GeV/cm'^ |50!], where the first uncertainty is from the 
slope of the circular velocity V{r) at the Sun's radius and the second is from uncertainty 
in distance of the Sun from the galactic center. The three studies are consistent with each 
other, but the first one finds significantly smaller uncertainty in the value of px- If we use 
the results of Ref. J48| , then the best fit regions plotted in Figs. [9] and [10] shift down by a 
factor of (0.4/0.3)^. Of course, it is not appropriate to simply scale up the signal by some px 
required to fit the positron and electron data; in a proper analysis, one should include the 
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appropriate likelihood for px and Vtxh"^ ■ We have also discussed some of the uncertainties 
in the velocity distribution function and their effects in Section lIVAt and, as with px and 
Vtxh"^-, these uncertainties should be included in the full likelihood analysis. 

Another effect that was pointed out recently |51[ is that the positrons produced within 
the unseen bound subhalos out in the Milky Way dark matter halo could propagate to 
the local neighborhood and contribute to the PAMELA and Fermi signal. For the best 
fit point, this study found that about 30% of the positrons in PAMELA could be due to 
annihilations in the subhalos |5l|. Taking px = 0.4 GeV/cm^ and an enhancement factor 
of 30% from subhalos, the discrepancy is reduced from a factor of 15 to a factor of 4. We 
see that these changes are insufficient to reach the best fit regime, even if all the parameters 
are simultaneously pushed in the most optimistic direction. A more troubling aspect of 
this calculation is however the fact that the subhalos considered were those resolved by the 
Via Lactea II simulation. This implies that the much larger contribution from the lower 
mass subhalos would over produce the signal. Appealing to subhalos to 0{1) changes to the 
required S'eff is therefore unmotivated. 

Finally, we may also appeal to small scale structure in the immediate local neighborhood 
to enhance the positron signals. Such a contribution is largely unconstrained by data, but 
present simulations, even without the deleterious effect of the disk of stars, do not predict 



large enhancements from such an effect [52[. A recent study [53[ using the Via Lactea II 
simulation and including the Sommerfeld effect for annihilation in the local clumps found 
that the probability of finding a subhalo close enough to make an 0{1) effect on the e^e~ 
flux above about 100 GeV is about 4% and that this contribution could be about 15% if the 
Via Lactea II subhalo mass function is extrapolated down to lower (unresolved) masses. Two 
points are worthy of further note here. First, these subhalo results are going to be crucially 
affected by the disk. The zeroth order expectation is that the gravity of the stellar disk will 
not allow such subhalos to survive or even if they do, truncate it signiflcantly. Second, for 
small m^ and near resonances self-interactions are important |2|, |4l| and that could modify 
the internal structure of the subhalo as well as the mass function. 

A more complicated set of astrophysical uncertainties is related to astrophysical back- 
grounds and the signals themselves. For example, the required Sommerfeld enhancement 
may be reduced if conventional astrophysical backgrounds have been under-estimated, cos- 
mic ray propagation models are modifled, pulsars contribute signiflcantly to the positron 
flux, or the PAMELA e^ sample includes some proton contamination. Any of these effects 
could reduce the required dark matter signal, but they could also plausibly ehminate the 



need for dark matter altogether, as has been argued in numerous studies |54|-l59l] . It is there- 
fore typically unclear what the motivation would be for assuming that such effects account 
for part of the signal but not all of it, independent of a desire to leave room for a dark 
matter signal. 

A possible exception to the above argument is the possibility of adjusting cosmic ray 
propagation models to match the Fermi spectrum. Such adjustments could remove the need 
to explain Fermi, without changing the PAMELA excess signiflcantly [60]. As shown in 
Fig. Uni for irix ^ TeV and m^ ^ 1 GeV, the allowed Sommerfeld enhancements may 
be in marginal agreement with PAMELA data. This scenario is, however, problematic. A 
contribution to the e^ flux at about 100 GeV to explain the observed PAMELA positron 
fraction implies a modiflcation to the e^ flux observed by Fermi. For example, a 500 GeV 
dark matter explanation for PAMELA with the correct relic density is marginally possible 
for m^ = 1 GeV according to the Ap flts of Ref. |2lJ as shown in Fig. [TOl This would. 
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however, imply a drop in the e^ flux beyond 500 GeV, which is not seen by Fermi, a point 



stressed in previous works [2l|, |46|. Analysis of such a scenario with modified propagation 



models would be another avenue for further study. 



C. Cosmic Microwave Background Bounds 



In this section, we include bounds from the effect of residual annihilation at last scattering 
on the ionized electron fraction and therefore on the cosmic microwave background (CMB) 
power spectrum. The CMB bound may be written as 6l|, l62 



.|„„.^ ^±^,}^(I^^, (35) 



e^ 



cos(./24/e 



/ VTeV 



where / is the average fraction of the energy produced in annihilation that reionizes Hydro- 
gen between redshifts of 800 and 1000. For the e^ final state, / ~ 0.7, while for the /i^ final 
state, / ~ 0.25 [62]. The left hand side S\v=o is the saturated Sommerfeld enhancement 
obtained by setting f = in Eq. ([6]). This is sufficient to approximate the Sommerfeld en- 
hancement during recombination, unless one is right on top of a resonance, which is highly 
disfavored by this bound. 

Equation (15BI) is not a monotonic function of ax- To impose this bound, we consider a 
region in ax that is bounded by values 10% larger and smaller than the value of ax dictated 
by the relic density bound of Qxh"^ = 0.114. We find that large resonances are no longer 
allowed by the CMB bound as a comparison of Fig. [10] and Fig. [11] will show. 



D. Force Carrier Decay Channels 

The best fit regions of the (mx, Scs) plane shown in Figs. |9] and [TO] assume that the force 
carrier decays solely to muons, leading to the XX -^ (jxp ^ Afi annihilation channel. The 4/x 
final state has been shown |20l . l2ll . l46l . |63| to be the final state that yields the best fit of all 



considered so far, as it leads to smooth e^ distributions that can simultaneously explain both 
PAMELA and Fermi. At the same time, this is an inefficient mode, as significant energy 
is lost to neutrinos. To ameliorate the discrepancy between the required Sce to explain the 
anomalies and the maximal S^s allowed by the relic density constraint, one might consider 
other channels, which could potentially allow similar contributions to the e^ spectrum with 
lower Ses and at lower mx- 

The existence of other decay channels is theoretically well-motivated and there are many 
possibilities. If the decay — )■ fi~^fi~ is possible, the decay — )■ e^e~ is also kinematically 
allowed. For m^ ~ 1 GeV, other decay modes, such as — )■ vrvr, KK, are also possible. As 
an example, if decays are controlled by its kinetic mixing with the standard model photon, 
then there will also be e^ and vr^ channels with roughly similar branching ratios, depending 
on m^ |6J]. Alternatively, particles may decay preferentially to the heaviest available 
states if they are coupled through standard model Yukawa couplings, or may have other 
interesting dynamics. In fact, in full generality, the particles may decay not only to pairs 
of standard model particles, but also through final states involving hidden Higgs bosons and 
other hidden particles, leading, for example, to XX annihilations to 6 or 8 particle final 
states. 
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FIG. 11: The maximal effective Sommerfeld enhancement Ses (black solid curve) after including 
bounds from CMB (see Sec. IV Cp compared to fits for ttt,,^ = 1 GeV for exclusive 4e and 4^ modes. 
The dark red (vertically hatched) region is the 99% contour for the 4/x mode and the light red 
(horizontally hatched) contour is for the 4e mode [2l(]. The best fit regions for the 4e mode are at 
lower 5efr and mx, but the best fit 4e mode point has a significantly lower likelihood than the best 
fit 4:fi mode point 2l| for a fixed cosmic ray propagation model. 



Studies so far have found that the 4/i mode is a better fit than other channels [2l|, |63 . 
A precise quantitative statement requires inclusion of correlations through covariance ma- 
trices for each experiment, and these are not available. Given this fundamental uncertainty, 
conclusions about the /x mode and other channels are far from firm and require additional 
assumptions. 

In Ref. 631], the 4e mode is shown to have a x^/dof that is larger by about ~ 1 compared 
the 4/i mode. These values of x^ are computed by assuming that the data points are 
uncorrelated and allowing for the spatial diffusion parameter to vary. Other results obtained 
by fixing the diffusion parameter seem to agree qualitatively with this result J2l|. In fact, 
marginalizing over rux and other parameters, Ref. 2l| found that the possibility of mixed 
modes in the ratio of 4e : 4/i = 1 : 1 is excluded at 99.9% CL. 

Despite these results, we may consider what impact other modes might have, especially 
in light of the fact that there is a lot of freedom in the cosmic ray propagation model 
parameters. As an illustration, we consider the best fit regions from Ref. J2l| for both the 
4e and 4/i modes separately. We plot these best fit regions with the relic density and CMB 
upper limits on Ses in Fig. [TTl The best-fit point for the 4e mode is at lower Ses and mx 
and the edge of the region is in marginal agreement with the relic density constraint. Of 
course, as discussed above, the best-fit point for the 4e mode has a much lower likelihood 
than the best-fit point for the 4/t mode, according to Ref. 2l|. For a 1:1 ratio, the best fit 
point would move to X masses somewhere between 1.5 and 2 TeV. This is disfavored by the 
relic density requirement even if we set the local dark matter density to be 0.4 GeV/cm^. 

For the reasons noted above, all of these conclusions are subject to many uncertainties. 
The present conclusion in the literature seems to that significant deviations from the pure 
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4/i channel do not provide a good fit to data. One way forward would be to redo fits of 
the kind discussed in literature, allowing for variations in cosmic ray propagation model 
parameters to see if this conclusion holds up in general. In any case, a careful statistical 
analysis is required before appeals to alternative channels alone may be considered as a way 
to allow Sommerfeld scenarios to explain the data. 

E. Bound States 

For the case of attractive forces, pairs of dark matter particles may form bound states 
by radiating force carrier particles if the force carrier mass m^ is below the binding energy 



a^mx/4 [l8|, |65|, |66|. As in the case of positronium, when this is kinematically possible, 
annihilation through bound states is preferred. In the limits m^ <^ Q;^r7ix/4 and v 3> 
iTi(f)/mx, the annihilation rate enhancement is about 20axv^{3 — f5)/(2f) for fermions, 
where v^ is the velocity of the final state particle when the bound state forms [l8|. This 
is roughly a factor of 7 larger than 5*" when f^ ~ 1. 

Unfortunately, such bound state effects do not help Sommerfeld enhancement models 
explain the PAMELA and Fermi data for two reasons. First, bound state formation is 
kinematically forbidden in large parts of parameter space that yield the best fits to data, 
and the large factor of 7 enhancement alluded to above is not realized for most of the 
parameter space plotted in Fig. ^U\ To see this, we first need to compute the value of ax 
required to get the right relic density. Now, bound states may form in the early Universe 
when they are energetically viable, that is when v < ax, and hence affect the relic density 
calculation. For simplicity we neglect this effect here, and use the same ax values from 
Fig. m which leads to a conservative bound. Then, for m^ = 0.25 GeV (1 GeV), the bound 
state enhancement formula quoted above is valid only for mx ^ 1-5 TeV (2.5 TeV). These 
inequalities are not satisfied for parameters that are most promising for fitting the data 
(m^ > 0.25 GeV and m^ < 4 TeV). 

Second, if bound states form, they are likely to be in confiict with bounds from the 
CMB. To see this, let's saturate the bound state enhancement at f ~ m^/irtx and esti- 
mate the enhancement at recombination as 20ax^™x/''^(/)- This should be compared to the 
CMB bound of (120//) (m^ /TeV) ~ 500 (m^/TeV), where we have set / = 0.25 to be 
conservative. To satisfy the CMB bound, then, we require ax < (m^/GeV)/40. This is 
generically violated if ax is fixed by the relic density. For example, for m^ = 1 GeV, the 
bound implies ax < 0.025. Requiring the correct relic density then implies mx < 1 TeV. 
However, we need mx > 2.5 TeV for bound states to form with m,^ = 1 GeV. The situation 
is no different with m^ = 0.25 GeV. Thus our saturation approximation above (for bound 
state enhancement) suggests that regions of parameter spaces where bound states can form 
are generically in conflict with the CMB. 



F. Non-minimal Particle Physics Models 

A potentially more promising direction is to construct more complicated particle physics 
models to achieve Sommerfeld enhancements that are large enough to explain PAMELA and 
Fermi. Such models will generically include additional annihilation channels. We begin by 
discussing the impact of these channels in general, and then examine various strategies one 
might explore to achieve larger Scs- 
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1. Additional Annihilation Channels 

As discussed above, to maximize Ses, we have considered only XX — )■ (pcf) annihilation, 
leading to the typical tree-level cross section estimate of Eq. ([2]). Even in the simplest 
models, however, one may have additional annihilation channels. For example, if is a U(l) 
gauge boson with mass generated by spontaneous symmetry breaking through the Higgs 
potential |A(|if p — w^)^, where if is a complex scalar, there is a physical hidden Higgs 
boson h with mass nih ~ \/2\v, where v is related to the 4> mass by m^ = y/Wnoxv. For 
perturbative A, m^ ^ 10 m^; h may be much lighter than 0, but it cannot be much heavier. 
The hidden Higgs boson therefore cannot be decoupled either dynamically or kinematically, 
and there is necessarily an additional annihilation channel XX — )■ </)/i. In the non-relativistic 
limit, and assuming there is no Yukawa coupling hXX, the cross section is 

a{XX -^ (f)h)v,ei = —^ ^^^^ 2W' ' (^6) 

fT^x ^^x (4m jf — m^^Y 

where 

I - I _ [(4m^ - (m^ + mhf) {Am\ - ("^0 - mnff'^ , . 



In the limit m^, rrih ^ fnx, 



a{XX ^ cph)vrei = 7^ . (38) 

4 rrix 



Additional annihilation channels and their impact on relic densities in Sommerfeld enhanced 



models have also been considered in, for example, Ref. [27 . 

To roughly accommodate such modifications, we may parameterize the effect of additional 
annihilation channels by assuming a tree-level annihilation cross section 



ira 



CTan^^reOo = k ^ , (39) 
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X 

where fc is a constant. If the new structures do not significantly modify the Sommerfeld 
enhancement S, the desired relic density will be achieved by modifying ax -^ ax/y/k, 
moving the positions of resonances and reducing the maximal 5'eff by the factor y/k. In 
general, Sommerfeld enhancement in the early Universe is important and the reduction in 
maximal Ses is smaller. The effects of varying k are shown in Fig. [121 

Of course, in general, new annihilation channels may have different couplings ax, differ- 
ent kinetic decoupling features, and different Tnt, all leading to different Sommerfeld effects, 
resonances in different places, and many other effects that cannot be captured by a single 
extra parameter k. The availability of additional annihilation channels also has a compli- 
cated effect on the signal. For example, for the case of decays to Higgs bosons and assuming 
h ^ 41 through real or virtual pairs, the h(f) mode leads to six or more leptons. This 
channel changes the energy spectrum of positrons that would be seen by PAMELA and 
Fermi and thus requires new fits to the data. The full analysis is therefore complicated, 
and requires a dedicated study incorporating all of these effects into the freeze out analysis. 
The main point is that non-minimal models will typically include additional hidden sector 
states and additional annihilation channels, generically reducing the maximum possible Sgg 
consistent with thermal relic density constraints. 
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FIG. 12: The maximal effective Sommerfeld enhancement factor S^s for k = 1,4, where the 
tree-level annihilation cross section is (can'WreOo = kira'x/m'x, for m^p = 0.25 GeV (left) and 
m0 = 1 GeV (right). The effects scale roughly as 1/Vk, with the deviation from 1/Vk scaling due 
to the Sommerfeld effect in the early Universe. 



2. Multi-state Dark Matter 



Two-state dark matter models have received considerably attention recently. For example, 



"inelastic dark matter" has been motivated by the DAMA Nal signal |67l . |68| and "exciting 
dark matter" has been proposed to explain the INTEGRAL/SPI excess towards the galactic 



center [69| . They have also been proposed as explanations of the PAMELA excess and the 

In these models, the stable dark matter 



WMAP haze towards the galactic center [38|, l70 



state is accompanied by a more massive unstable state, which is separated from the stable 
state by a mass gap that is chosen to fit the data and is much smaller than the dark matter 
mass. The presence of this almost degenerate excited state significantly complicates the 
early Universe freeze out analysis. The Sommerfeld enhancements are also different |25;] . To 
compute the relic density, one would need to include the additional annihilation channels 
and the two-state dark matter Sommerfeld enhancement in a self-consistent early Universe 
freeze out calculation. 

Dark matter models with multiple stable dark matter states have also been pro posed to 



explain the PAMELA and Fermi excesses and the DAMA/LIBRA signal [7ll-l73| . These 
models provide an interesting avenue to reduce the required Sommerfeld enhancement to 
explain both the PAMELA and Fermi excesses 73j. This is because the heavier (~ TeV) 
component contributes mainly to the Fermi excess, whereas the lighter (~ 100 GeV) dark 
matter particle contributes predominantly to the PAMELA positron excess. The required 
Sommerfeld enhancements are about 100 for the TeV state. 

Such enhancements may appear to be consistent with our relic density constraint to within 
a factor of 2, but the relic density calculation that we have presented here is not complete 



in these scenarios. In the simplest such models [73| , the two states couple to the dark gauge 



boson with the same charge. As stated, however, these models are ruled out by the direct 
detection limits, because the Xq — )• Xq cross section due to the kinetic mixing would be 



too large for e > 10 ^ [17|, |36|. The way out of this constraint is to specify that each state 



is further split into a metastable and stable state with very small mass gap and that the 
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particles only have off-diagonal couplings, so that kinetic energy larger than the mass gap is 
required to initiate Xq — )■ X*q, where X* denotes the excited state of X jlTJ]. These mass 
gaps may then be chosen to be ~ C(MeV) or ~ C(0.1 MeV) to explain INTEGRAL/SPI and 
DAMA/LIBRA signals respectively [73|. The relic density calculation is now significantly 
more complicated. Additional annihilation channels, for example, annihilations between the 
heavy and light states and annihilation to h(j), and the presence of an almost degenerate state 
must be taken into account in the early universe calculation. The Sommerfeld effect for the 
lighter state is generically important because the coupling to the dark gauge boson would be 
set by the relic density constraint for the heavy particle, and so ax ~ 0.02, much larger than 
the required ax for thermal dark matter with mass ~ 100 GeV. The new channels, such 
as the h(j) channel, also modify comparisons to PAMELA and Fermi as these new channels 
may contribute to the e^ flux with a different energy spectrum. These issues are beyond 
the scope of the present work. 

Multi-state dark matter theories may also have significantly suppressed k. For example, 
in the non-Abelian models of Ref. [7^, values of k as small as 0.14 are possible, and these 
have recently been claimed to provide marginally consistent explanations of the cosmic ray 



anomalies [75| . In these studies, however, the irreducible annihilations to Higgs bosons have 
not been included, and the Sommerfeld effect on freeze out has been neglected. Given the 
low value of k, these changes to the annihilation cross section may be important. It would 
be interesting to see if the consistency remains after a complete analysis. 

To summarize, finding a model with more degrees of freedom that requires smaller boosts 
to fit the PAMELA and Fermi data does not by itself imply an improvement, as the bounds 
on S'eff will generically also be stronger. Generalizations of the freeze out analysis described 
here are required to determine if these scenarios may self-consistently explain PAMELA and 
Fermi data. 



3. Very Heavy Force Carriers with m^ ^ 1 GeF 

One reason the maximal S'eg is somewhat limited is that the resonances are not significant 
for m^ ^ 1 GeV and m^ ^ 1 TeV. For larger m^, however, the resonances move to larger 
mx, as evident from Eq. ([7]). Generically, for m^ > 2 GeV, decay may produce anti- 
protons. However, if such contributions to the anti-proton spectrum may be otherwise 
suppressed or accommodated [76|, such large m^ may allow one to improve the consistency 
of Sommerfeld-based explanations of the PAMELA and Fermi data. 

As an example, in Ref. [l9(], a leptophilic model was constructed where the anti-proton flux 
is suppressed by dynamics, freeing the force carrier mass to be much larger than ~ GeV. 
For 777,0 ~ 10 GeV and 777^ ~ 800 GeV, the authors noted that resonant Sommerfeld 
enhancement could enhance signals with enhancement factors as large as 5* ~ 1000. However, 



the study of Ref. [19j did not include the Higgs annihilation channel, and did not include the 
impact of Sommerfeld enhancement on freeze out. Note also that, as discussed in Sec. IIIIBl 
for very large 777^ near the freeze out temperature, particles may not be in complete thermal 
equilibrium with the standard model, modifying the freeze out analysis. Nevertheless, such 
possibilities certainly merit further study. Note, though, that these models are inaccessible 
to searches for GeV forces. 
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4- Running ax 

In all of our analysis above, we have assumed that the running of the gauge coupling 
in the dark sector is not significant. Although this is certainly a good approximation for 
the minimal Abelian model we have considered, the situation is much more complicated 
for non-Abelian models. It has been recently pointed out that a non-Abelian hidden sector 
could explain both PAMELA and Fermi, while resulting in the correct relic density |77| . 
This work notes that the Sommerfeld effect is a soft process and the enhancement factor 
S is determined by the coupling at scale m^, whereas for the tree-level cross section the 
coupling must be evaluated at nix- For a non-Abelian theory where the coupling constant is 
smaller at higher energies, this could significantly relax the tension between the relic density 
constraint and the requirements for explaining PAMELA and Fermi. 

However, this analysis |77| did not include the Sommerfeld effect on the annihilation in the 
early Universe. This would be particularly important because the Sommerfeld enhancement 
is now much larger even in the early Universe. The model also includes a Higgs doublet 
H and the dark matter annihilation XX — )• HZ must also be included. The stability of 
the Higgs sector is a major concern for these models because the Higgs annihilation cross 
section is too small to suppress their relic density |77|. The decay lifetimes that have been 
computed [77| are long (~ 10^ s or longer) and could destroy the successful predictions of 



Big Bang nucleosynthesis (BBN) because of the injection of electromagnetic energy [78|, |79 . 
Given these issues, it is not clear that these models can successfully satisfy the relic density 
constraints, but future work addressing these points could be an interesting way forward to 
explain PAMELA and Fermi excesses using Sommerfeld-enhanced annihilation. 

We note that the arguments regarding the decay of the hidden sector Higgs boson are 
quite general. The hidden Higgs boson should decay to standard model particles, otherwise 
it would contribute to the dark matter abundance. In doing so, the decays should not 
destroy the successful predictions of BBN. 

VI. CONCLUSIONS 

In this work, we have investigated scenarios with Sommerfeld-enhanced annihilation of 
thermal dark matter and determined the largest possible indirect signals, subject to the most 
basic of constraints, namely, that the thermal dark matter have the correct thermal relic 
density. In addition to their potential impact on cosmic e^ spectra, such scenarios may also 
predict observable effects in the cosmic microwave background and the diffuse extragalactic 
7-ray background, excess 7-ray, neutrino and radio wave emission from the galactic center, 
and excess 7-rays from dwarf spheroidal satellites of the Milky Way and extended regions 



centered on the galactic center J2l|, |46|, |6l|, l62l. |8Ch88| . The absence of each of these potential 



signatures provides significant constraints. The bounds derived here are complementary to 
those, in that these bounds are derived by studying not the annihilation products, but the 
production. These results are therefore model-independent and robust, in that they assume 
only that the dark matter is produced by thermal freeze out, which is the central assumption 
motivating the consideration of Sommerfeld enhancements in the first place. 

This analysis included the possibility of resonant Sommerfeld enhancement and the effect 
of Sommerfeld enhancement in the thermal relic density calculation. The possibility of 
Sommerfeld enhancement at low velocities, with annihilation cross sections proportional to 
v~J, or even v~J on resonance, implies that the thermal relic density is sensitive to the dark 
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matter density's evolution long after conventional freeze out at Tf ~ rnxf^b. As a result, 
the relic density depends on many aspects absent in conventional scenarios. These include 
the cutoff of resonant Sommerfeld annihilation, the equilibration of force carrier particles, 
the temperature of kinetic decoupling T^d, and the efficiency of self-scattering to thermalize 
the dark matter velocity distribution after T^^- 

These dependences may be very important. A non-generic but intriguing example is that, 
on resonance, dark matter freezes out but may then later melt in, chemically recoupling for 
a second era of efficient annihilation. As a result, exact resonances suppress rather than 
enhance indirect signals, in contrast to naive expectations. More generally, the Sommer- 
feld enhancement of freeze out implies that smaller tree-level annihilation cross sections are 
required to give the correct relic density, and so reduces the effective Sommerfeld enhance- 
ment of indirect signals. We have found reduction factors of a few, with several ~ 10 — 30% 
variations possible, depending on the value of T^d and other freeze out parameters. 

For the minimal scenario analyzed here, and adopting the most optimistic assumptions 
to maximize indirect signals, the largest possible values of Ses are 7, 30, and 90 for m^ = 
250 MeV and mx = 0.1, 0.3 and 1 TeV, respectively. As shown in Figs. 1^ and [TU| such 
enhancements fall short of explaining the PAMELA and Fermi cosmic e^ excesses. For the 
best fit point, as discussed in Sec. |Vl this discrepancy is over an order of magnitude and 
cannot be eliminated by appeals to enhanced local dark matter density or boosts from small 
scale structure. Bound state effects are also unable to eliminate the inconsistency, as they 
are typically insignificant in the parameter regions that may potentially explain the cosmic 
ray anomalies, and are in any case stringently constrained by the CMB. Along with the 
astrophysical uncertainties discussed in Sec. |V]and the complementary constraints from the 
cosmic microwave background and other sources, these results imply that the dark matter 
motivations for Sommerfeld enhancements and GeV dark forces at present are, at best, 
strained. 

We have discussed several possibly interesting directions for further study in Sec. |Vl 
Non-minimal models will generically include additional annihilation channels, such as those 
involving hidden Higgs bosons. These will generically strengthen the bounds on Ses- We 
noted also that the decay of the hidden sector Higgs boson to standard model particles 
could be constrained by demanding consistency with the measured light element abundances 
produced during BEN. Additional features and constraints particular to specific models have 
been noted in Sec. El We have focused on the 4/i final state here, which generically give the 
best fits to the data. Different annihilation channels, such as 4e final states, are optimized at 
lower Ses and mx, but also yield worse fits to the data according to present estimates in the 
literature. However, this conclusion may depend significantly on the cosmic ray propagation 
model parameters used. At present, no Sommerfeld models have been shown to explain 
the PAMELA and Fermi excesses consistent with the constraint of thermal freeze out, and 
we stress that complete and dedicated studies of freeze out are required to judge properly 
whether more complicated models can self-consistently explain PAMELA and Fermi data. 

The maximal Sommerfeld enhancements derived here, although not sufficient to explain 
current PAMELA and Fermi excesses, are nevertheless significant enhancements. These will 
be probed as PAMELA and Fermi continue to gather data, and, in the near future, at AMS 
on the International Space Station. The interpretation of future data as dark matter may, of 
course, continue to be clouded by confusion from astrophysical uncertainties, but the diverse 
set of particle physics and astroparticle physics experiments that may probe weak-scale dark 
matter in the near future at least provides some reason for optimism. 
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